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We would like to determine \Vcb\ from the exclusive semi-leptonic decay B D*lv. The 
differential decay rate is 

^ = ^{w^ - lY'^m'j,,{mB - mn*fg{w)\V,b\^\rB^D,{w)\\ (1) 

where w = v ■ v' and ^(1) = 1. At zero recoil {w = 1) heavy-quark symmetry requires 
-?^B^D*(1) to be close to 1. So, \Vcb\ is determined by dividing measurements of dV/dw by the 
phase space and well-known factors, and extrapolating to w 1. This yields \Vcb\J^B^D*{^)-, 
and J-B^D*{^) is taken from "theory." To date models [ ^ or a combination of a rigorous 
inequality plus judgment [ |2| have been used to estimate Tb^d*{^) — 1- In this work [ ^] we 
calculate J^b-^d*{^) with lattice gauge theory, in the so-called quenched approximation, but 
the uncertainty from quenching is included in the error budget. 

The "form factor" J^b^d* {w) is a linear combination of several form factors of the matrix 
elements (L'*|V^|-B) and {D*\A^\B). At zero recoil all form factors but are suppressed by 
phase space, so 

Tb^D'{1) = /iAi(l) = {D*{v)\A^\B{v)), (2) 

which should be "straightforward" to calculate in lattice QCD. But a brute force calculation of 
{D*\A^\B) would not be interesting: similar matrix elements like {Q\A^\B) and I(k\V^\B) have 
15-20% errors [|, |]. 

Thus, we have to involve heavy-quark symmetry from the outset: if we can focus on hA-^ — 1, 
we have a chance of success, because a 20% error on — 1 is interesting: 0.2 x 0.1 = 0.02. 
There are three specihc obstacles to overcome: (i) statistical uncertainties, (ii) normalization 
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uncertainties in the lattice axial vector current, and (iii) how to treat heavy quarks since 
m{,a 5^ 1. The first two need computational insight; the last two theoretical insight. In the last 
several years, we have developed tools to attack these problems [ ^, |7[ |8|, |9|, [Tl| . 
At zero recoil heavy-quark symmetry implies [ |l^ 



^Ai (1) — riA llsgur-Wiso + OLukc + + ^l/m^ 



(3) 



where r/^ is a short-distance coefficient of HQET, and the di/m" contain long-distance matrix 
elements. The structure of the ^/toq corrections is 
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where S = l/(2mc) -|- l/(2m;,) and A = l/(2mc) — l/(2mft). One must make sure to calculate 
r]A and the ^s in the same renormalization scheme. 

Lattice gauge theory with Wilson fermions has the same heavy-quark symmetries as con- 
tinuum QCD, for all mqa. It therefore admits a description with HQET, provided mq ^ A [ 
^, ^ |9|, [l^. In this description, HQET matrix elements, such as the is in Eqs. (^) and (|5|), 
are essentially the same as for continuum QCD. So, one needs some quantities with small sta- 
tistical and normalization errors, whose heavy-quark expansion contains the ^s. Then, one 
calculates the short-distance part in perturbation theory [ 14], extracts the £s from a fit, and 
reconstitutes hAiiX)- 

In our work on the B ^ D form factor [ 11], we found that certain ratios have the desired 
low level of uncertainty: 



{D\c-i%\B){B\lri^c\D) 
{D\c'y^c\D){B\b-f%\B) 
{D*\c-/'^b\B*){B*\h'^c\D*) 
{D*\cj^c\D*){B*\b-f%\B*) 
{D*\c-i^-l^b\B){B*\h-ii-^r^c\D) 
{D*\cj^j5c\D){B*\h^j5b\B) 



{??1.-*[i-^pA2-43)a2s]}', 
l-iAA^-i^^^A^y . 
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For lattice gauge theory, the heavy-quark expansions in Eqs. (P)-(^) have been derived in Ref. [ 
0], leaning heavily on Refs. [ ITH]. The one-loop expansions of r^y* and ?7]|* are in Ref. J p!0{ . 

(3) 

Thus, these ratios yield all terms in ^i/^s except i\j . 

We wish to obtain the l/wig corrections to the double ratios, but the lattice action and 
currents do not normalize all such terms correctly. HQET reveals several sources of such 
contributions, in a systematic way [ 15, P]. The most crucial are the I/tuq corrections to the 



currents, which enter the double ratios as follows: 



[1 - \{Xh/ml - l/mcrrib + Xe/mg)]^ 
[1 - A(2A, - l)/m2][l - A(2Xfe - l)/ml] 
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where A is proportional to Ai or A2, and Xq/mQ indicates incorrect normalization. The correct 
normalization of the l/rricmij terms is built into the current we used. The cancellation of the 
others is a key feature of the double ratios. Other contributions either vanish or are correctly 



2 



normalized to order [ ^. This, and other matching uncertainties of order and (A/mg)'^, 
are put into the error budget. 

Fig. ||(a) shows the heavy-quark mass dependence from Ref. [ |3| . As expected, iy is the 
largest of the I/^tt-q matrix elements. Because of the fit, the value of iy is highly correlated 

(3) 

with that of iy , but the physical combination is better determined. 

We have also studied the dependence of the calculation on the mass of the light spectator 
quark, over the range 0.4 < mq/nis < 1. As seen in Fig. 1(b), there is a slight linear dependence 
on m^, which is proportional to nig. The points are correlated, so the trend is significant. The 
main effect of extrapolating in is to increase the statistical error. In addition, there must 
be a pion loop contribution [ !£], which is mistreated in the quenched approximation. We treat 



the omission of this effect as a systematic error. 
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Figure 1: (a) Heavy-quark mass dependence of the double ratios, (b) Chiral extrapolation of 
/^Ai(l). 



3 



After putting everything back together again, we find [ ^] 



T n ^ — n Qi q+0.024 _|_ n ni c+0.003 +0.000+0.006 /I 

J-B^D* - U.yid.Q 017 ± L).Uib_o oi4-0.016-0.014> (J-Uj 

where the uncertainties stem, respectively, from statistics and fitting, HQET matching, lattice 
spacing dependence, the chiral extrapolation, and the effect of the quenched approximation. In 
Fig. ^(a) we compare our result for Tb^d* (1) against the estimate based on the quark model [ 
1^], and on the sum rule The defects are as follows: The quark model omits some dynamics 
(more than quenching), and it is not clear that it gives the £s in the same scheme as rjA- The 
sum rule has an incalculable contribution from excitations with (M — m/j*)^ < /x^, which can 
only be estimated. The present lattice result is in quenched approximation, but the error from 
quenching is the last error bar in Eq. (p!o[). 

For using this result in a global fit to the CKM matrix, it is useful to have some idea what 
values of theoretical quantities are more (or less) likely. A flat distribution would be incorrect, 
because the first error in Eq. ( p!o| ) is essentially statistical, and the others are under some control. 
Also, one cannot rule out a tail for lower values; they are just unlikely. Finally, we know that 
JFb^£)*(1) < 1 [|2|. A simple formula that captures these features is the Poisson distribution 

P{x) = Nx'e-'\ ^^ ^-:Fb^d*{1) _ ^^^^ 
^ ' ' 0.087 ^ ' 



In the future one could reduce the uncertainty by a factor of 3, as sketched in Fig. |2|(b), and 
one could provide a distribution stemming from the Monte Carlo calculation, and properly 
propagated through the systematic analysis. 
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Figure 2: (a) Comparison of methods, (b) Simple Ansatz for more (and less) likely values, now 
and in the future. 
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